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Abstract: Landscape structure is often regarded as an important factor 
that governs the distribution and abundance of species. Therefore it is 
critical to understand the landscapes and their dynamics. Patterns of 
landscape elements strongly influence the ecological characteristics. This 
study was designed to document and map the current status of the tropi¬ 
cal dry deciduous forest of the Tadoba-Andhari Tiger Reserve (TATR), 
Central India, (using IRS P6 LISS IV data) and to describe its landscape 
structure at three levels of organization viz. landscape, class, and patch. 
The study area was classified into 10 land cover classes that include 6 
vegetation classes. The landscape structure was analyzed using FRAG- 
STATS using 12 set of indices. The TATR landscapes have a total of 
2,307 patches with a mean patch size of 25.67 ha and patch density of 1.7 
patches per km 2 . Amongst all land cover classes, mixed bamboo forest is 
dominant—it occupied maximum area (77.99%)—while riparian forest 
is least represented (0.32%). Mixed forest has maximum number of 
patches among all vegetation classes. Results have shown that despite 
being dominant in the area, mixed bamboo forest has low patch density 
(0.25/100 ha). Dominance of mixed bamboo forest is attributed to large 
patch sizes and not to the number of patches. This study has focussed on 
the approach of integrating satellite forest classification and forest inven¬ 
tory data for studying forest landscape patterns. 
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Introduction 

The management of wildlife and protected areas has several 
goals: conservation, optimal use of forest resources, meeting the 
demands for the scientific underpinnings of managing landscapes, 
and incorporating the consequences of spatial heterogeneity. The 
sensitivity of ecological effects of resource management towards 
spatial configuration is gaining acceptance worldwide. The land 
cover classification from remote sensing data is a powerful tool 
that can provide repetitive and spatial information concerning the 
landscape (Chust et al. 2004). 

Studies have highlighted the role of remote sensing data from 
earth observation satellites in forest ecology (Jadhav et al. 1990; 
Innes and Koch 1998; Skole and Tucker 1993 and Franklin et al. 
1994). Satellite remote sensing plays a crucial role in generating 
information about forest cover, vegetation types and land use 
changes (Cherill and McClean 995). Both coarse and fine spatial 
resolution satellites have proved to a significant resource for 
forest mapping. Milanaova et al. (1999) studied broad vegetation 
type stratification using coarse resolution data like NOAA- 
AVHRR and several studies have also reported mapping at finer 
resolution data of LANDSAT TM+ (Groom et al. 1996; Guillem 
et al. 2004). Remote sensing estimates of regional variation in 
biodiversity can be used in analyzing diversity patterns, land¬ 
scape patterns and aiding conservation efforts (Gould 2000). 
Spatial patterns have a strong influence on the information con¬ 
tent of ecosystem components. Since landscape structure is often 
regarded as an important pre-requisite that governs the distribu¬ 
tion and abundance of species, the first step is to understand 
landscapes and their dynamics. It is widely acknowledged that 
patterns of landscape elements strongly influence the ecological 
characteristics. Therefore, spatial pattern characterization and 
quantification of land-cover classes to relate pattern and process 
is a pre-requisite at landscape level (Turner 1987). It is not only 
important to understand how much there is of a particular com- 
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ponent but also how it is arranged (Turner 2001). Quantification 
of landscape pattern is necessary for understanding the composi¬ 
tion and configuration of landscapes. The underlying premise is 
that the explicit composition and spatial form of landscape mo¬ 
saic affect ecological systems in ways that it would be different, 
if the mosaic composition or arrangement were different (Wiens 
1995). Moreover, recent studies have demonstrated that land use 
and landscape changes significantly affect biodiversity (Cousins 
and Eriksson 2002; Gachet et al. 2007; Miyamoto and Sano 
2008). Spatial tools of remote sensing and GIS have a capacity to 
quantify land-cover patterns and understand spatial heterogeneity 
(Turner 1990; Sinha and Sharma 2006). Methods are needed to 
quantify aspect of spatial patterns that can be correlated with 
ecological processes (O’Neill et al. 1988). 

A large number of spatial indices are based on patch metrics 
that quantify the spatial pattern at three different levels of or¬ 
ganization—patch, land cover, and landscape—using the pro¬ 
gramme FRAGSTATS (McGarigal and Marks 1995). Numerous 
studies have advocated the authenticity of this spatial pattern 
analysis programme (Lu et al. 2003; Cushman et al. 2008; Lele 
et al. 2008). The landscape structure of the U.S. state Kansas was 
analysed at three scales by calculating the landscape pattern met¬ 
rics (Griffith et al. 2000). The utility of landscape pattern indices 
was studied for judging the habitat implications of alternate land¬ 
scape plans or designs (Correy et al. 2005). LISS III data of IRS 
ID and IRS P6 were studied and compared for landscape eco¬ 
logical application by computing class level and landscape level 
metrics (Sharma 2009). 

Previous studies conducted in the study area (Dubey & Mathur 
1999; Dubey 1999) were restricted to the classification of vege¬ 
tation using IRS LISS II data of 36.5m coarse resolution. How¬ 
ever, neither of the studies have ever dealt with exploring the 
landscape structure of TATR. Satellites with high spatial resolu¬ 
tion facilitates detailed assessment of vegetation, identification of 
smaller patches and is also considered helpful in evaluation of 
impacts on biodiversity of specific management policies (Innes 
and Koch 1998). IRSP6 LISS IV is one such high resolution 
satellite. Several studies have assessed the utility of LISS IV in 
different thematic areas: Gupta and Jain (2005) carried out urban 
mapping; Kumar and Martha (2004) assessed damage from land¬ 
slide; Rao and Narendra (2006) mapped and evaluated urban 
sprawling; Kulkarni et al. (2007) studied glacial retreat in Hima¬ 
laya; and Bahaguna (2004), Sarangi et al. 2004; Rajankar et al. 
(2004) evaluated application of IRS P6 in coastal and marine 
zones. Since this satellite it is expected to provide a new oppor- 
tunitiesy to efficiently make detailed vegetation cover maps effi¬ 
ciently for large study areas, the present study was initiated with 
the aim to document and map the current status of the TATR 
forest. The IRS P6 LISS IV data was of 5.8 m high resolution 
and was used to describe the landscape structure of TATR at 
three spatial levels of organization: viz. landscape level, class 
level, and patch level. The primary objective of thise study is to 
understand the habitat quality of TATR for better conservation 
management of the tiger reserve for large mammals by providing 
general characterization of the landscape structure and patterns. 


Materials and methods 

Study area 

The Tadoba-Andhari Tiger Reserve (TATR) represents an im¬ 
portant habitat for wildlife in Maharshtra State in Central India. 
It comprises of Tadoba National Park (TNP) and Andhari Wild¬ 
life Sanctuary (AWLS).The area lies between 20°04' to 28°025'N 
and 79°13' to 79°33'E (Fig. 1). It covers 625.40 km 2 and has a 
subtropical climate with three distinct seasons: summer, mon¬ 
soon, and winter. The climate is characterized by a hot and pro¬ 
longed summer while winter is short and mild. TATR is mostly 
undulating and hilly in the north and almost a plain in the south¬ 
ern part of the TATR. It has been classified as southern tropical 
dry deciduous forest (Champion and Seth 2005). 

Data source 

Three digital scenes of IRS-P6 LISS IV with 5.8m resolution 
were acquired from the National Remote Sensing Agency 
(NRSA) for December 2005 and January 2006. During this pe¬ 
riod the vegetation was in full bloom and cloud free data could 
be obtained. 

Data analysis 

Image preparation for classification 

Unwanted artefacts like additive effects due to atmospheric scat¬ 
tering were removed through set of pre-processing or cleaning up 
routines. First-order radiometric corrections were applied, using 
the dark pixel subtraction technique (Lilesand & Kiefer 1994). 
This technique assumes that there is a high probability that at 
least a few pixels within an image should be black (0% reflec¬ 
tance). Because of atmospheric scattering, the imaging system 
records a non-zero Digital Number (DN) value at the supposedly 
dark- shadowed pixel location. Therefore the DN value was sub¬ 
tracted from the data to remove the first-order scattering compo¬ 
nent. Images were then registered geometrically. Uniformly dis¬ 
tributed Ground Control Points (GCPs) were marked with a root 
mean square error of less than one pixel and the image was re¬ 
sampled by the nearest neighbourhood method. All the scenes 
were mosaiced and the study area was extracted using digital 
boundaries. A reconnaissance survey was conducted from Febru- 
ary-June 2006 to have gain an overview of broad vegetation 
types in the study area and to devise a sampling strategy for de¬ 
tailed sampling analysis. An approach of systematic stratified 
sampling was used and range maps were used developed to strat¬ 
ify the area for ground-truthing and vegetation plots. Later, from 
June 2006 to February 2008, intensive ground- truthing was con¬ 
ducted and a total of 810 GPS points, including 520 vegetation 
plots, were collected to capture the variation in spectral signa¬ 
tures of different vegetation types over the entire study area and 
to achieve higher accuracy of vegetation mapping. Some 810 
spectral signatures were collected out of which 284 were used for 
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the accuracy assessment. The data collected from the plots was tion communities existing in the field, 
subjected to hierarchical cluster analysis to determine the vegeta- 


Tadoba-Andhari Tiger Reserve (TATR), Central India 



Fig. 1. Study area 



Image classification 

Unsupervised classification was done on image, using the near¬ 
est-neighbour algorithm for differentiating spectral reflectance of 
various objects (Lilesand and Kiefer 1994). It examines the simi¬ 
lar pixels in an image and aggregates them into a number of 
classes. Initially, the entire area was classified into 40 classes, 
which were iterated 10 times with convergence threshold of 0.98. 
According to Champion and Seth (2005), the area is classified 
under group 5 and subgroup 5A as Southern Tropical Dry De¬ 
ciduous Forest. The entire area was divided into 40 classes ini¬ 
tially. Later, these 40 classes were classified under the broad er 
classes of forest, waterbody, scrub, open forest and agricul¬ 
ture/settlement. The non-forest classes were masked and cluster 
analysis was conducted on the data acquired from vegetation 
plots to determine the forest communities. Results of cluster 
analysis were then used to classify forest classes. Finally, 40 
classes were merged into 10 classes, including six vegetation 
classes which were well observed on the ground. The pixellated 
classified output image was obtained. Since the objective of the 
study is to study habitat quality of TATR for large ungulates, 
therefore, map was subjected to two 3><3 filters to ensure 
smoothening and the patches below 0.5ha were removed. Area of 
0.5 ha was considered as insignificant to study habitat of large 
ungulates. Finally, the area was calculated for each class. The 
accuracy of the map was tested using some of the ground truth 


points not used during classification. The land cover information 
of these locations was compared to classified maps. Accuracy 
was tested using Kappa Statistics (Khat coefficient) (Lillesand 
and Kiefer 1994). 

Image preparation for Fragstats 

For Fragstats input, rasters 0s (background data) was recoded as 
negative integers so that it would be ignored by Fragstats. 

Landscape structure analysis 

For the quantification of the TATR landscape of TATR, uses of 
statistical measures or indices were made that describe landscape 
configuration and composition. These indices were calculated by 
FRAGSTATS (Mcgarigal and Marks 1995). The FRAGSTATS 
is a spatial pattern analysis programme for categorical maps. It 
simply quantifies the areal extent and spatial configuration of 
patches within landscape. There are two versions of FRAG¬ 
STATS, Vector (ARC/INFO) and Raster (Image Maps) versions. 
The raster version has been used in this study to compute metrics. 
The landscape structure was analyzed at three different spatial 
scales, viz. , landscape, class and patch levels, using 12 sets of 
indices as shown in Table 1. The reason to selected these indices 
were attributed to their simplicity. These were the basic indices, 
which provided the general description of the landscape. Numer¬ 
ous studies have supported their authenticity of these indices 
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(Griffith et al. 2000 and Cushman and Neel 2008). Fig. 2. de- scribes the methodology in the form of a flow chart. 


Table 1. Metrics used for the landscape characterization of TATR 


Level 

Metrics 

Description 

Unit 

LI Landscape 

No. of Patches (NP) 

No. of patches in a landscape 

None 

L2 Landscape 

Patch Density (PD) 

No. of patches in a landscape divided by total landscape area. 

No./lOOha 

L3 Landscape 

Largest Patch Index (LPI) 

Area of the largest patch in the landscape divided by total landscape area 

% 

L4 Landscape 

Interspersion and Juxtaposition (LJI) 

Adjacency among patches of different class 

% 

L5 Landscape 

Simpson Diversity Index (SIDD 

Diversity of patches in the landscape 

None 

L6 Landscape 

Simpson Evenness Index (SIEI) 

Even distribution of area among patch types 

None 

Cl Class 

Percentage of Landscape (PLAND) 

Percentage of landscape comprised of corresponding class. 

% 

C2 Class 

No. of Patches (NP) 

No. of patches of corresponding class. 

None 

C3 Class 

Patch Density (PD) 

No. of patches of corresponding class. 

No./lOOha 

C4 Class 

Mean Patch size (MPS) 

Average patch size of the corresponding class 

ha 

C5 Class 

Interspersion and Juxtaposition (LJI) 

Adjacency among patches of corresponding class 

% 

PI Patch 

Patch area (Area) 

AREA equals the area (m 2 ) of the patch divided by 10,000 (to convert to ha) 

ha 


Fig. 2. Flowchart of methodology 


Results 


Landscape composition 



Seven communities were identified from cluster analysis: Tec- 
tona, Tectona-Cliestanthus, Tectona-Chloroxylon, Zizyphus- 
Adina, Pterocarpus-Flacourtia, Dalbergia-Mitrigyna and Ripar- 

<£] Springer 


ian community. Mapping was done using the cluster analysis as a 
premise (Fig. 3). As a result, 10 landcover/landuse classes were 
delineated: teak forest, teak mixed bamboo forest, mixed forest, 
mixed bamboo forest, riparian forest, grassland, scrub, open 
forest, agriculture/settlements and water body (Fig. 4). Amongst 
the vegetation classes, mixed bamboo forest occupied the maxi¬ 
mum proportion of the study area (75.81%) and the riparian for¬ 
est occupied the least (0.61%). Table 2 shows the areas of the 
different classes mapped. 


Tcctonn - Diospyros 
Tcctona- Clicstanthus 
Tcctona- Chloroxylon 

Zizyphus-Adina 
Pterocarpus-Flacourtia 
Dalbcrgia-Mitragyna 
Syzigi um -1 lardwickia 



Fig. 3. Dendrogram showing the different communities using hierarchi¬ 
cal cluster analysis 

Table 2. Area of landcover classes delineated 


S.No. 

Class 

Area (ha) 

% Area 

1 . 

Teak Forest 

1385.4 

2.22 

2. 

Teak Mixed Bamboo Forest 

789.5 

1.26 

3. 

Mixed Forest 

3996.2 

6.39 

4. 

Mixed Bamboo Forest 

47402.4 

75.81 

5. 

Riparian Forest 

380.3 

0.61 

6. 

Grassland 

2717.4 

4.35 

7. 

Open Forest 

1585.9 

2.54 

8. 

Scrub 

3174.7 

5.08 

9. 

Agriculture/ S ettlement 

755.4 

1.21 

10. 

Water Body 

337.7 

0.54 


Total 

62525.4 

100 
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Landuse/landcover map of Tadoba-Andhari Tiger Reserve 


Fig. 4. Landuse/landcover map of Tadoba-Andhari Tiger Reserve 


Landscape pattern and structure 

Landscapes in the study area had been defined at three levels of 
hierarchy starting from broader levels to narrower ones, i.e. land¬ 
scape level, class level and patch level. 

At Landscape Level 

Tadoba-Andhari Tiger Reserve’s landscape was found to be 
heterogeneous in nature. As shown in Table 3, a total of 2307 
patches of different types and varying patch sizes (mean patch 
size =25.67 ha) could be recognized in the landscape with a 
patch density of 1.7 patch per km 2 . The landscape was evenly 
interspersed with different forest types as indicated by the inter- 
spersion value of 50. The landscape was not very diverse in na¬ 
ture as shown by low values of Simpson Diversity index 0.38 
and Simpson Evenness index 0.42. 


Table 3. Landscape metrics for TATR landscape 


Landscape Metrics 

Values 

No. of Patches (NP) 

2307 

Patch Density (PD) 

1.7/km 2 

Largest Patch Index (LPI) 

32.53% 

Interspersion and Juxtaposition (TJI) 

50 

Simpson Diversity Index (SIDI) 

0.38 

Simpson Evenness Index (SIEI) 

0.42 


At Class Level 

The landscape was composed of six major vegetation types, viz., 
mixed bamboo forest, mixed forest, teak forest, teak mixed bam¬ 
boo forest, riparian forest and grassland have described. The 
result of all the metrics computed at the class level is given in 
Table 4. 


Table 4. Class level metrics for landscape of TATR 


Vegetation Types 

PFAND 

(%) 

NP 

PD 

(No./1 OOha) 

MPS 

(ha) 

FPI 

(%) 

IJI 

(%) 

Mixed Bamboo Forest 

77.9 

340 

0.25 

136.1 

32.5 

68.2 

Mixed Forest 

6 

671 

0.49 

5.3 

0.6 

5.8 

Teak forest 

2 

182 

0.13 

6.6 

0.6 

61.6 

Teak Mixed Bamboo Forest 

1 

42 

0.03 

13.7 

0.2 

14 

Riparian Forest 

0.3 

35 

0.03 

2.3 

0.02 

62.8 

Grassland 

4.1 

225 

0.16 

7.2 

0.6 

42.8 


Area and density metrics: Amongst all vegetation types, 
maximum percentage of land was covered by mixed bamboo 
forest (77.99%) with the highest mean patch size (136.09 ha) and 
the largest patch index (32.53%). Mixed forest has the highest 
number of patches (671) and therefore has the highest patch den¬ 
sity (0.49/1 OOha). However, riparian forest had the least land¬ 
scape area (0.3%) with lowest number of patches (35) and lowest 
mean patch size (2.25 ha). The average patch size varies from 
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2.25 ha to 136 ha. Apart from mixed bamboo forest, all forest 
types have mean patch size below 15 ha. 

Interspersion Metrics: The value of the interspersion / juxtapo¬ 
sition metric was found to be highest in mixed bamboo forest 
(75.23%) followed by riparian forest (62.79) and teak forest 
(61.63%). The least interspersion was found in mixed forest 
(5.78%). 

At Patch level 

The area of each patch comprising a landscape mosaic is the 
most useful piece of information contained in a landscape. The 
analyses revealed that among all the vegetation classes, mixed 
forest has the maximum number of patches among all vegetation 
classes. The small-size patches ranging from 0.5 to 5 ha were 
highest in number (599) and the large-size patches (above lOOha) 
were few. Among all vegetation classes in the TATR landscape, 
riparian forest had least number of patches (35) with 18 patches 
of size ranging from 0.5-5 ha and only one big patch of 22 ha. 

Discussion 

Vegetation Mapping 

Landscape elements, coupled with satellite imagery can be effec¬ 
tively used to monitor biodiversity (Nagendra and Gadgil 1999). 
At high spatial resolution, many factors affect the recorded re¬ 
flectance of the plant communities (species, crown closure, 
crown geometry, stand density, soil moisture and sun angle). 
This made it possible to map the communities. Despite using 
high-resolution satellite data, problems were encountered while 
mapping land-cover classes, since the presence of bamboo in the 
understory created a spectral overlap with spectral signatures of 
crown cover. Special consideration was given to the compatibil¬ 
ity of collected ground data and the spectral qualities measured 
by satellite. Among the land-cover classes, mixed bamboo forest 
was the most dominant class in TATR because of extensive 
bamboo growth. Despite, teak being a dominant tree species in 
the study area, little area was occupied by pure teak patches. It is 
present only in the northern part of Tadoba National Park. Teak 
was mostly present along with its associates and bamboo in other 
forest types. The old plantations of teak have now been con¬ 
verted into teak mixed bamboo forest, because of extensive 
flowering of bamboo in TATR in the mid 1980s (Dubey and 
Mathur 1999). The scrub and open forest were mostly found in 
the southern part of the Andhari Wildlife Sanctuary. It could be 
attributed to villages present in the periphery of the southern 
zone, which exert an anthropogenic pressure that leads to the 
degradation of the surrounding forests. The presence of natural 
water sources and high protection status are major reasons con¬ 
tributing to the presence of all vegetation classes in Tadoba Na¬ 
tional Park. 

Landscape Characterization 

Structural analysis of the landscape helps in problem identifica¬ 
tion and the degree of severity, which is useful in ecosystem 


management (Formon and Godron 1986). The analysis here sup¬ 
ports the observation that a small set of indices can capture sig¬ 
nificant aspects of landscape pattern. Structural analysis of the 
TATR landscape reveals its heterogeneous nature with large 
variations in patch size, but with low diversity, low evenness, 
and intermediate interspersion of forest types. Mixed bamboo 
forest covers the maximum area of TATR and therefore, the 
maximum percentage of the TATR landscape (77.99%), indicat¬ 
ing its dominance among vegetation classes. Percentage land¬ 
scape at the class level is a good indicator of fragmentation 
(Sinha and Sharma 2006). Mixed forest was found to be the most 
patchy, with the highest number of patches (671) and the highest 
patch density (0.49/100 ha). Our results have shown an interest¬ 
ing pattern, however: despite being dominant in the area, mixed 
bamboo forest has low patch density (0.25/100 ha). Despite hav¬ 
ing few patches, it had highest mean patch size (136 ha). Domi¬ 
nance of mixed bamboo forest is attributed to large-size patches, 
not to the number of patches. 

Mixed bamboo forest and riparian forest had the highest adja¬ 
cencies among all the vegetation types, indicating that these two 
forest types share their edges with the rest of the forest types. 
Nevertheless, teak mixed bamboo forest and mixed forest had the 
least interspersion among all forest types, due to their clumped 
distribution in the landscape. 

This study has focussed on the approach of integrating satellite 
forest classification and forest inventory data for studying forest 
landscape patterns. IRS P6 LISS IV data have proved to have 
immense potential for minutely capturing the structural details of 
the landscape due to its high resolution and multispectral nature. 
This attribute has also been used for analyzing the patch dynam¬ 
ics in the landscape. The results presented here support focusing 
on a few metrics that represent overall landscape structure for 
landscape characterization and monitoring. Park managers 
should consider indices as tools for comparing different land¬ 
scapes patterns. The trends depicted by the application of land¬ 
scape metrics may be assimilated into their prognostic models 
and scenarios to support strategic decision making for regional 
conservation and policy development. 

Conclusion 

The TATR landscape is comprised of 10 landuse/landcover 
classes. Six major vegetation types—mixed bamboo forest, 
mixed forest, teak forest, teak mixed bamboo forest, riparian 
forest and grassland—were delineated. Mixed bamboo forest was 
the most dominant vegetation class, covering 75.81% of TATR, 
and riparian forest was the least represented, covering 0.61%. 
Detailed landscape-level analysis on the number of patches, 
patch characteristics (composition), spatial arrangement, and 
proximity of different patches (configuration) has provided cru¬ 
cial information on the landscape structure. Structural analysis of 
the TATR landscape revealed its heterogeneous nature with large 
variations in patch size. The landscape was found to have uneven 
distribution of the patches with intermediate interspersion among 
forest types. The dominance of mixed bamboo forest is attributed 
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to the large size of the patches, despite being less in number. 
Mixed forest was found to have highest number of patches (671) 
and therefore the highest patch density (0.49/ha), while riparian 
forest has lowest patch density (0.03/ha). The results indicate 
that the landscape metrics in FRAGSTATS are effective in char¬ 
acterizing the landscape. This component of the our study has 
demonstrated the efficacy of high-resolution satellite IRS P6 
LISS IV data with multispectral capability in detailed mapping 
of forest types with sharp boundaries and accurate area estimates. 
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